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Abstract 

We show that grand unified theories based on SO (10) generate quite naturally 
baryon number violating dimension seven operators that violate (B — L), and lead 
to novel nucleon decay modes such as n — > e~K + , e~7r + and p — > vir + . We find that 
in two-step breaking schemes of non-supersymmetric SO( 10), the partial lifetimes 
for these modes can be within reach of experiments. The interactions responsible for 
these decay modes also provide a new way to understand the origin of matter in the 
universe via the decays of GUT scale scalar bosons of SO(10). Their (B — L)-violating 
nature guarantees that the GUT scale induced baryon asymmetry is not washed out 
by the electroweak sphaleron interactions. In minimal SO (10) models this asymmetry 
is closely tied to the masses of quarks, leptons and the neutrinos. 
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Baryon number violation is a very sensitive probe of physics beyond the Standard Model 
(SM). Interactions which violate baryon number (B) are not present in the renormalizable 
part of the SM Lagrangian, but can arise through effective higher dimensional operators. 
The leading B violating operators [1] have dimension six and are hence suppressed by two 
powers of an inverse mass scale. These operators arise naturally when SM is embedded 
in grand unified theories (GUTs) such as SU(5) and SO (10). They lead to nucleon decay 
modes such as p — > e + n° and p — > VK + , which conserve baryon number minus lepton 
number (B — L) symmetry. Experimental searches to date have primarily focussed on 
these modes with the latest limits on proton lifetime constraining the masses of the heavy 
mediators to be larger than about 10 15 GeV. This is in accord with the scale inferred from 
the unification of gauge couplings. 

Going beyond the d = 6 baryon number violating operators, the next-to-leading ones 
have d = 7, and obey the selection rule A(B — L) = —2 for nucleon decay [2]. These 
operators lead to novel nucleon decay modes such as n — > e~K + ,e~7r + , and p — > vk + , 
which have received less attention. In this Letter we show that these d = 7 operators arise 
naturally in unified theories based on 50(10), upon the spontaneous breaking of (B — L), 
which is part of the gauge symmetry. In particular, we find that in non-supersymmetric 
SO (10) models with an intermediate scale so that gauge couplings unify, the partial lifetime 
to these decay modes can be within reach of ongoing and proposed experiments. Further- 
more, we show that these new modes provides a novel way to understand the origin of 
matter in the universe. This mechanism relies on the fact that, owing to their (B — L) 
breaking nature, a GUT scale induced baryon asymmetry would not be affected by the 
electroweak sphalerons [3] and would survive down to low temperatures. Observed baryon 
number of the universe then carries the direct imprint of GUT scale physics. This is unlike 
the (B — L)-preserving baryon asymmetry induced in the decays of GUT mass particles 
such as in SU(5), which is however washed out by the sphaleron interactions, leaving no 
trace of GUT physics. We show that in minimal 5*0 (10) models [4] which have been highly 
successful in predicting large neutrino oscillation angles, including a relatively large value 
of sin 2 2# 13 ~ (0.085 — 0.095), consistent with recent results [5], the baryon asymmetry of 
the right magnitude is generated by the new (B — L)-violating mechanism. The results of 
this paper should provide motivations to search for (B — L)-violating semi-leptonic decay 
modes of the nucleon in the ongoing and the next round of searches. Their observation 
would furnish evidence against the simple one-step breaking of GUT symmetry, and could 
also resolve the mystery behind the origin of matter in the universe. 

We start by writing down the d = 7 5-violating effective operators in the SM [2] in the 
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standard notation for fermion fields: 
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Here stands for the covariant derivative, and if (1,2, 1/2) is the Higgs doublet. These 
operators obey (B — L) = +2 selection rule and mediate nucleon decays of the type n — > 
e~K + ,e~7r + and p — > un + . We first show how these operators arise naturally in 5*0(10) 
theories [6] when (B — L) symmetry contained in it is broken. This breaking may occur at 
the GUT scale as in models with supersymmetry, or at an intermediate scale Mj below the 
GUT scale, as in non-supersymmetrci SO(10), which requires such a scale to be compatible 
with gauge coupling unification. To see the origin of Eq. ([I]) in SO(10) via scalar boson 
exchange, we write down the Yukawa couplings in the most general setup. Noting that the 
fermion bilinears contain 16 • 16 = 10 s + 120 a + 126 s , the Yukawa couplings are [8]: 



£(1646,10*,) = hi 



[u c i Q j + v'Lj) h - (d^Qj + e ?Lj) h + i^-QiQj + u\e] - d 



+ [euld] + Q l L j ) co c ] 
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with h and / being symmetric and g being anti-symmetric in flavor indices The 
SU(3)c x SU(2)l x U(1)y quantum numbers of the various sub-multiplets in Eqs. (j2D- 
& are: h(l, 2, +1/2), h(l, 2, -1/2), w(3, 1, -1/3), w c (3, 1, 1/3), p(3, 2, 1/6), p(3, 2, -1/6), 



3 



V (3, 1, 2/3), r7(3, 1, -2/3), $(3, 3, -1/3), $(3, 3, 1/3), X (3, 2, 7/6), *(3, 2, -7/6), 5(3, 1, -4/3), 
and 5(3,1,4/3). 

The (B — L) generator of 50(10) is broken by the VEV of the SM singlet field A c 
in 126#, which has (B — L) = —2. This VEV supplies large Majorana masses for the 
right-handed neutrinos through the coupling fij\Z6i/fUjA c . It also generates trilinear scalar 
couplings of the type p*uH, r]*pH, p*$if and x*vH, thereby inducing the d = 7 baryon 
number violating operators of Eq. (pQ), via the Yukawa couplings of Eqs. ([2])-(jl]). The flavor 
symmetric Yukawa couplings of Eqs. ©-© generate the operators 3 and 0\ through the 
diagrams shown in Fig. [TJ 



Q L u c L 




(AO H (AO H 

(a) (b) 

Figure 1: Effective baryon number violating d = 7 operators induced by the symmetric 
Yukawa couplings of 10# and 126# of 50(10). 

The trilinear couplings in Fig. [Tjhave different sources in 50(10). The quartic coupling 
(126) 4 , which is invariant, contains the term (2, 2, 15) • (2, 2, 15) • (1, 1, 6) • (1, 3, 10) under 
SU(2) L x SU(2) R x SU(A) C subgroup. The p*{3~, 2, -1/6) C (2,2,15), #(1,2,1/2) C 
(2, 2, 15), while u(3, 1, —1/3) C (1, 1, 6), and this coupling would contain the term p*ojHA c . 
The three non-trivial invariants of the type (126) 2 • (126*) 2 also contain this trilinear term. 
Similarly, the three quartic couplings (120) 2 • (126) 2 would generate the trilinear terms 
p*7]H*, x*vH, and §*pH* vertices, and along with Eq. (JJJ) would induce the remaining 
non-derivative d = 7 operators of Eq. ([TJ. For a detailed discussion see Ref. [5]. 

Such trilinear couplings are also present when the 126// is replaced by a 16# albeit 
in a slightly different way. The 16 h contains a SM singlet filed with B — L = +1 which 
acquires a GUT scale VEV, and a h(l, 2, —1/2) field with B — L = —1. The trilinear scalar 
couplings of the type 16#16/fl0# and 16#16#10h would mix the B — L = Higgs doublet 
h(l, 2, 1/2) from the 10 # and the h(l, 2, 1/2) Higgs from the 16 H which has B — L = +1. 
The light SM Higgs doublet then would have no definite B — L quantum number. The 
(1, 2, 4) component of 16# contains the field p*(3, 2, —1/6), and the (2, 1, 4) of 16# contains 
u(3, 1, —1/3), and thus the coupling p*uH is generated via the 16#16#10# coupling. 

The d = 7 operators of Eq. flTJ can also arise by integrating out the vector gauge boson 
Vq(3,2, 1/6) and V^((3, 1, -2/3) of 5O(10), which lie outside of SU(5) [6j. The covariant 
derivative for the 126# would contain the term VqV u cH(A c )^ which generates the d = 7 
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operators. When 16# is used instead of the 126#, the covariant derivative would contain a 
similar term, but now with (B — L) = +1 and —1 for H and A c C IQh respectively. 

Partial lifetime for (B — L) violating nucleon decay: The diagrams of Fig. [1] lead to the 
following estimate for n — > e~n + lifetime^] 

r(B e -^)«dH„ ^^(1 + D + Ff^ {^fj^-y (5) 

Here we have defined the Yukawa couplings of u and p fields appearing in Fig. [T] to be 
Yqq w and YhdFp- The factors D and F are chiral Lagrangian factors, D ~ 0.8 and F ~ 0.47. 
f3 h — 0.012 GeV 3 is the nucleon decay matrix element, Vr = (A c ), and v = (H°) = 174 
GeV. We have defined the trilinear coupling of Fig. [TJto have a coefficient Xvr. The mass 
of u)(3, 1, —1/3) is constrained to be relatively large, as it mediates d = 6 nucleon decay. 
For Y w 10 -3 , M w > 10 11 GeV must be met from the d = 6 decays. As an illustration, 
choosing Y QQu = Y Ld c p = lO" 3 , M u = 10 11 GeV, M p = 10 8 GeV, and Xv R = 10 11 GeV 
in Eq. ()5]), we find r„Ri 3 x 10 33 yrs. Such a spectrum is motivated by the intermediate 
symmetry SU{2) L x SU{2) R x SU(4) C , which is found to be realized at M 7 w 10 11 GeV 
from gauge coupling unification [7]. As a second example, take M p = 10 6 GeV, M w = 10 16 
GeV, Xvr = 10 16 GeV, Yqq^ = Y Ld c p = 3 x 10~ 3 . This choice of spectrum leads to 
T n ~ 4 x 10 33 yrs. This spectrum can arise as follows. Suppose the p(3, 2,1/6) particle, 
along with a pair of (1, 3, 0) scalar particles (contained in the 45#, 54# or 210# needed for 
symmetry breaking) survive down to M/ = 10 6 GeV. The SM gauge couplings are found 
to unify at a scale Mx ~ 10 15 GeV in this shown in Fig. |2j This scenario would 

predict observable rates for both the (B — L)-conserving and (B — L)-violating nucleon 
decay modes. Analogous results are obtained from the exchange of 77(3, 1, 2/3) — p(3, 2, 1/6) 
scalar bosons from 120 h- Since these particles do not induce d = 6 baryon number violation, 
they can both have mass of order Mj, which would enhance the nucleon decay rate. 

While the AB = —2 nucleon lifetime is in the experimentally accessible range for 
reasonable choice of parameters as shown, it is quite sensitive to the precise values of the 
intermediate scalar masses. For example, a factor of 3 increase in M p and M w will increase 
the lifetime by a factor of 10 4 . Not finding these modes will not exclude this class of 5O(10) 
models, but a discovery of the AB = —2 nucleon decay mode would lend strong support 
to a new mechanism of baryogenesis via (B — L)-violating decays of scalars, to which we 
now turn. 

New baryogenesis scenario at the GUT epoch: We now present a new baryogenesis 
scenario at the GUT epoch, using the (B — L)-violating decay of the scalar uj (3, 1, —1/3) 
with a GUT scale mass. The magnitude of the asymmetry is directly linked to the neutrino 
masses, since the Yukawa couplings that induce the asymmetry are the same couplings 



1 (-B+L)-preserving nucleon decay has been studied in the context of i?-parity breaking SUSY in Ref. [9]. 
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Figure 2: Unification of the three SM gauge couplings obtained with a light p(3, 2, 1/6) and 
two (1,3,0) scalar multiplets at Mj = 10 6 GeV. 

that are involved in neutrino mass generation. (B — L asymmetry in decays of specific 
heavy particles has recently been discussed in Ref. [ID].) To be concrete, we shall work in 
the framework of non-supersymmetric 5*0(10), although our results would hold for SUSY 
S'O(IO) as well, with some minor modifications. The couplings of Eq. d2])-(H]) imply that 
u has two-body decays into fermions of the type u —> QQ, u c e c , v c d c , u c d c , Q L. These 
decays preserve (B — L), as can be seen by assigning (B — L){oj) = —2/3. Now, u also has 
a two-body scalar decay, u — > pH* as shown in Fig. [3] (a), which uses the (B — L) breaking 
VEV of A c . The scalar field p has two-body fermionic decays of the type p — )■ Ld c , v c Q 
(the latter if kinematically allowed), which define [B — L) charge of p to be +4/3. Thus 
the decay u — > pH* would violate (B — L) by —2. 




(c) (d) 

Figure 3: Tree-level diagram and one-loop corrections responsible for generating (B — L) 
asymmetry in u decay. 
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Let the branching ratio for u — » pH* be r which produces a net (B — L) number of 
4/3, and that for — > p*H be r, with net (B — L) = —4/3. The branching ratio for 
the two-fermion decays u — > ff is then (1 — r) which has (B — L) = —2/3, and that for 
(V* — / / is (1 — r) which has (B — L) = 2/3. Thus in the decay of a u> + u* pair, a net 
(.B — L) number, defined as €b-l, is induced, with €b-l = (B — L) w + (B — L) w * = 2(r — r) 
The loop diagrams for u> —¥ pH* are shown in Fig. |3] (b)-(d), which involve the exchange 
of fermions. Since u can also decay to two on-shell fermions, these loop diagrams have 
absorptive parts and also CP violation. 

We evaluate Fig. [3] in a basis where the Majorana mass matrix M v c of the v c fields is 
diagonal and real. The contributions of Fig. [3] (b)-(d) to ^b-l are found to be 

Br r f 

t B -L = - , A , 2 Im \\v R Tr |yJ Xa;S Y Qu cp M v . F{M U , M p , M v o) Y V . LH 
-Y d o u c w Y\ cLp Y u c LH M u c F(M p , M u , M v] ) 

-YL^> Y ^ F '( M ^ M -'' M i)( X ' v RT}} > (6) 

where the three terms are in order from Fig. [3] (b), (c) and (d). Here we have defined the 
trilinear scalar vertices of Fig. |3] (a) and (d) to have a coefficients Xvr and X'vr in the 
Lagrangian. Yq v <^ is the Yukawa coupling matrix corresponding to the coupling Q z/ c p, etc. 
Br stands for the branching ratio Br(u; — > pH*). A factor of 2 has been included here for 
the two SU(2) l final states in the decay. The functions F and F' are defined as 

F(a,6,c) = ln(l + a 2 /c 2 ) + 9(l-c 2 /6 2 )(l-c 2 /6 2 ), 

F'(a,b,c) = (l-c 2 /a 2 )/(l-6 2 /a 2 )0(l-c 2 /a 2 )(l-c 2 /a 2 ). (7) 

Here O stands for the step function, signalling additional ways of cutting the diagram when 
Mj < M p or Mj < Mu in Fig. |3l Fig. [3] (d) arises because in any realistic 50(10) model 
there are at least two u fields. The heavier u field is denoted as u'. We have also assumed 
that M u — M^ ^> T w , so that there is no resonant enhancement for the decay. 

To estimate Br = Br(u; — > pH*) appearing in Eqs. (Q, let us assume that u is the field 
u from 10# with Yukawa couplings as given in Eq. 02]). The partial widths for the decays 
Yi(u — > pH*) and T 2 (p —> f f) are then given by 



(8) 



with Br = r 1 /(r 1 + T 2 ). For M w = 10 16 GeV, h 33 = 0.6 (corresponding to the top quark 
Yukawa coupling at GUT scale) with other negligible, and \vr = (10 14 , 10 15 , 10 16 ) GeV, 
one gets Br = (1.4 x 10" 4 , 1.4 x 10" 2 , 0.58). 

The (B — L) asymmetry es-L of Eq. (J6]) will result in a baryon to entropy ratio Yb 
given by 

Y B = = 6 -^d , (9) 

s g* 
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where g* = 130 is the total number of relativistic degrees of freedom at the epoch when 
these decays occur. The factor d in Eq. ([9]) is the dilution factor which takes into account 
back reactions that would partially wash out the induced baryon asymmetry. Defining 
K = r ^ u ~?£ H ^ , where H is the Hubble expansion rate, H = 1.66 gl ■ff-, the dilution 
factor can be written as p.] d ~ 1 (K < 1), and d ~ K ^ )a . e (K > 1). For M w = 10 15 GeV, 
and Xv R = (10 14 , 10 15 , 10 16 ) GeV, we find K = (0.12, 12.3, 1230) and the corresponding 
dilution factors to be d = (1.0, 1.4 x 10" 2 , 7.5 x 10~ 5 ), with Br = (1.3 x 10~ 2 , 0.58, 1.0). 

We now show how the GUT scale induced baryon asymmetry in u> — > pH* decay can 
consistently explain the observed value of Yb = (8.75 ± 0.23) x 10 -11 , in a class of minimal 
SO (10) models. In these models a single 10h and a single 126# couple to fermions, as in 
Eqs. (|2])-(j3]). It has been shown that these models lead to large mixing angles for solar and 
atmospheric neutrino oscillations. Furthermore, they predict sin 2 2^13 w (0.085 — 0.095), 
both in the non-SUSY and the SUSY versions [I], which is consistent with recent results 
from Daya Bay and other experiments [5]. To illustrate how realistic choice of parameters 
generate acceptable Yb, we choose the u field to be almost entirely in the 10#. We also 
choose X'vr that appears in Fig. 131(d) to be small, so that the leading contribution to cb-l 
is from Fig. El (c). In this limit, we find e B - L ~ ^-^^^{l +\n(l + Mf/M^)} sin0. Here 
we have kept only the third family Yukawa couplings, and defined <p — argl/if^/lA + f }■ 
Choosing h 33 ~ 0.6 (the top quark Yukawa coupling at the GUT scale), and A = 0.25, 
vr = 10 16 GeV, /3 = 10~ 2 (so that f 3 VR = 10 14 GeV, consistent with the light v T mass 
arising via the seesaw mechanism), <ft = 0.12, we find €r-l = 1-6 x 1.9 x 10 -5 . If M w = 10 15 
GeV, then Br = 0.96, K = 197 so that the dilution factor is d = 5.6 x 10~ 4 . This results in a 
net Yb = 8.2 x 10 -11 , consistent with observations. We emphasize the intimate connection 
between €b-l and neutrino masses, since Yq Lui ,, Yg u cp, etc present in es-L are the 126 
couplings that determine the neutrino masses via the seesaw mechanism. 

In conclusion, we have shown that all d = 7 baryon number violating operators that 
lead to nucleon decay modes such as n — > e~K + , e~ir + and p — > utt + , emerge naturally as 
effective low energy operators in a wide class of S'O(IO) models. In non-supersymmetric 
5*0(10) models with an intermediate scale, we find the rates for these nucleon decay modes 
to be within reach of experiments. We have also shown that the existence of these (B — L)- 
violating interactions allows a new scenario for baryogenesis where a (B — L) asymmetry is 
generated in the decay of GUT mass particles which survives to low temperatures unaffected 
by the sphaleron interactions. In minimal 50(10) models which predict large neutrino mix- 
ing angles, including dis, this new mechanism can explain the observed baryon asymmetry 
of the universe. 
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Foundation Grant Number PHY-0968854. 
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